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1. Introduction 

The current observational data Q strongly support Dark Energy domination in the Uni- 
verse. The nature of the Dark Energy is still a mystery and the "coincidence question" 
unanswered. So, alternative pictures exist in the literature, among which is the one pro- 
posed recently in the framework of Brane- World cosmology ||, |||, ||, based on the obser- 
vation that due to energy exchange with the bulk, the present accelerating Universe may 
be a late-time stable attractor of the cosmological evolution equations || . This is an effort 
towards a natural resolution of the "coincidence problem" of cosmology, which in addition 
leads to several phenomenologically interesting properties at the fixed point. However, to 
accommodate the correct amount of matter Q m ~ 0.3 at the fixed point, one has to assume 
that the brane has negative tension, with the exponential expansion on the brane driven 
by energy influx from the bulk ]?]]. 

However, negative tension branes are believed a priori to be unsatisfactory for a realistic 
physical model, for essentially two reasons. First, it is known that the gravitational force 
on a negative tension brane in Randall-Sundrum-like scenarios [|| (with AdS bulk) becomes 
repulsive. Indeed, following H one can show that the four-dimensional Newton's constant 
Gn in this case is proportional to the brane tension a with positive definite coefficient, 
leading to negative Gn for negative a. Additional matter in the bulk, for instance in the 
form of a minimally coupled scalar field with normal kinetic term and arbitrary potential, 
does not improve the situation ||. Second, negative tension RS-II branes are believed to 
be unstable. 

However, the physics changes considerably, if one extends the model to include the Ricci 
scalar term on the brane, induced by matter quantum loops or by finite brane-thickness 
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effects. The action of this, so called DGP, model || in its simplest version (without bulk 
cosmological constant or matter) is 

S =\ j d i xdy v ^M 3 R + J d 4 Xy/=g~M 3 AK + ~ J d A x^/=gM 2 R + a J d^x-s/^g (1.1) 

and describes a 3-brane with tension a but no extra matter on it, embedded in Minkowski 
4+1-dimensional bulk. R is the aforementioned intrinsic scalar curvature term on the 
brane, while AK = K + — K~ is the jump of the trace of the extrinsic curvature across the 
brane. As usual, the purpose of this modified K dependent Gibbons-Hawking (GH) term 
is to cancel unwanted terms in the variation of the bulk action related to the discontinuities 
of the derivative of the metric across the brane, so that one ends up with the proper Israel 
matching conditions 1 . 

The solution for the metric, relevant to cosmology which is our main interest here, is 
found to be 

/ „, ,x9 ( -1 



Sab 



V o i 

where Z2 symmetry across the brane is manifest and e = ±1. The parameter H is given 
by 

# = 3t * E+ ' =* H = ^L±Jl + °J£-\. (1.3) 
M 2 3M 2 M 2 I y 3M 6 J 

Metrics with opposite H are related by time reversal. So, only two of these four values of 
H correspond to independent solutions. We choose the ones which correspond to the two 
values of e = dbl, both with the + sign in front of the square root. 

We are interested in solutions with the following properties: (a) H > 0, in order to 
correctly describe the present accelerating expansion of the Universe, (b) eH < in order 
to improve the chances for stability, (c) attractive Newton's law on the brane, and (d) all 
the above consistent even with a < 0. This, may eventually allow for a natural explanation 
of the cosmic acceleration and a resolution of the "coincidence issue" H . 



Notice, though, that none of the solutions given in (1.2), ( |1.3j ) satisfies these require- 
ments. Closest to being satisfactory is the self-accelerating one, corresponding to e = +1. 
This has H > even with a not too negative a, but the metric in this case grows in the bulk 
and leads, perhaps not surprisingly, to ghost instabilities in the spectrum of perturbations 
1 11, [l^, [13|, [l4|, |l^, pill . For a positive tension brane a helicity-0 excitation of the spin-2 



graviton is a ghost, while for a negative tension brane the spin-0 mode becomes a ghost. 
For tensionless brane the ghost field is a linear combination of the spin-0 mode and the 
helicity-0 excitation of the graviton. 

The other independent solution with e = —1 is stable, but it is unsatisfactory in 
connection with (d) above, because one needs positive cosmological constant on the brane 
to explain the accelerating cosmic expansion. Is there a way to satisfy all four requirements 
(a)-(d)? 



1 A very explicit derivation of the variation of the GH term can be found in juj. 
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A step in this direction was outlined in O 2 . The proposal was to modify the bulk 



action (1.1) by the multiplication with an appropriate smearing function near the brane as 
follows 

Sbulk = \J d 4 xdyy/^§M 3 R - ~ J d 4 xdy^M 3 F(x,y)R (1.4) 

where in the simplest case F(x,y) = 1 — m5(y) with m a parameter and 8{y) a 5-like 
function with a second parameter a. A simple choice is 5{y) = TT~ 1 a/(a 2 + y 2 ). The 
extrinsic curvature terms should be modified accordingly. The model still has the solution 
Ql.2| ) for the metric, with H now given by 

'-g^-di'^ 1 + afl.(i- S )' J- <L5) 

The solution with the + sign in front of the square root, e = — 1 and in the limit a —* 
and m — > with m/(ira) ~ const > 1, has H > 0, describes a self-accelerated brane, which 
satisfies all our requirements, with the exception of the stability issue which cannot be 
decided at the level of such a non-covariant formulation of the model. Nevertheless, as it 
was shown in such a modification leads to a flip of the sign in front of the extrinsic 
curvature terms in the brane equations of motion. Exactly this sign is responsible for the 
stability of the metric perturbations and the absence of unstable modes in the normal 
branch. Thus, modifying the action in the above mentioned way we expect to obtain self- 
acceleration while keeping equations as they are in the normal stable branch. Further nice 
properties of this modification are considered the relaxation of the bulk gravity scale and 
the conservation of the number of parameters. 

The purpose of the present paper is to develop a manifestly covariant formulation of 
the above modification, making use of a bulk scalar field instead of the function F{x,y). 
Naturally, extra complications arise from the fact that one has to satisfy also the scalar 
field equation of motion, as well as the corresponding additional matching condition on the 
brane. 



2. Codimension-1 brane in 5-dimensional bulk 

Consider a 3-brane embedded in a 5-dimensional bulk. The coordinates in the bulk are 
denoted by x A with capital latin indices running from to 4 and with x 4 = y. The 
coordinates on the brane are denoted by with Greek indices from the middle of the 
alphabet taking values from to 3. Occasionally, we shall use t instead of the coordinate 
with index 0, while spatial indices on the brane will be denoted by lowercase Latin letters 
a, b, ... . The position of the brane in the bulk is parameterized as x A = X A (^), X 4 = Y. 
Finally, we shall be using tildes to designate quantities referring to the bulk. Thus, the 
bulk metric is c/ab and the induced metric on the brane is = d tl X A d v X B gAB{X A (£ i fJ ')). 
Here = d/d^ and our convention for the signature is (— , +, +, +, +). 

2 For a related approach see also JL7J. 
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We shall be interested in the model described by the generic action 

Sbulk = J d 4 xdy^ (M1f(<I>)R - ^X{$)~g AB d A <&d B $ - , (2.1) 

5branc = / d*Xy/=g (^M 3 F($)AK + L hiane ( 9flv , VO) . (2.2) 

The bulk scalar field <J? is supposed to be a modulus field. This is the field relevant to the 
hereby proposed modification of the DGP model. V"(3>) may have a constant term A, being 
the cosmological constant in the bulk. F, X and V do not depend on g. AK = K + — K~ 
is, as in the jump of the trace of the extrinsic curvature across the brane. Standard 

Model fields ip confined on the brane may also be present. They will not be needed in our 
discussion. We choose to study the equations of motion (EOM) in the Gauss-Normal (GN) 
coordinate system, with the brane located at y = 0. In GN coordinates it is convenient to 
identify X 11 = so that and are indistinguishable (static gauge). Further, in this 
coordinate system g = diag((^, 1) and the induced metr ic on the brane is g^ v — 9fiv\y=o- 
The non- vanishing Christoffel symbols are = — ^^f^ 1 , = \(j vp ^%^- and f%/, where 
r%/|y=o = Fftu, the being the Christoffel symbols obtained from the brane metric g^ u . With 
x^ and ^ identified, we have = djdx^ = d/d^, while d A = d/dx A . 

As a side remark, notice that if F(&) is not a constant, it can be brought to any 
other non-constant form by an appropriate redefinition of the field The two actions are 
equivalent at the classical level. A similar observation applies to the function X(<&). If it 
is non- zero, it may be transformed to any non- vanishing constant. The two redefinitions 
cannot, in general, be applied together. 

2.1 Bulk 



Varying the bulk action (2T) with respect to the metric one obtains in the bulk 



- M 3 FG AB + M 3 (D A D B F - ~g AB UF) + Xd A $d B $ - ^~g AB ~g CD d c <S>d D <S> - g AB V = 0. 

(2.3) 

Here G AB = R AB — \Rg AB is the bulk Einstein tensor, D is the covariant derivative for 
the bulk metric g AB and □ = g AB D A D B . Similarly, the EOM of the scalar field <3? reads 

XO§ + ±—g AB d A <f>d B $ + - = 0. (2.4) 

Here the superscript ^ denotes the n-th derivative with respect to the field 

For simplicity we shall restrict ourselves to diagonal bulk metrics. Furthermore, fol- 
lowing the cosmological principle, we shall take space on the brane to be homogeneous and 
isotropic. In the special case of zero spatial curvature on the brane the most general ansatz 
for the metric in GN coordinates is then 

/ -N 2 (t,y) \ 

~9AB= A 2 (t,y)5 ab . (2.5) 

V o o i ; 
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The non-zero Christoffel symbols for the (/x, v) part are t\ t = N/N, = AA5 a b/N 2 , 
T at = A5 b a /A. In what follows, we will use dot for time derivative and prime for y derivative. 



Direct calculation of the Einstein tensor gives the following non-zero components 



G t 



G a b 



(-iV 2 )3 
{A 2 



A 2 



A' 



A 2 N 2 + A 2 + 

J_ (P A 
' N 2 \A 2+ A 



A^_ 
~A 



) An\ AT_ A" n A'N' N^_ 
'AN + ~A? + ~A + AiV + ~N 



G 



G, 



ty 



2.2 Brane 




(2.6) 
(2.7) 
(2.8) 
(2.9) 



Life on the brane is described by the equations of motion on the brane. They are obtained 
from the brane lagrangian, supplemented by the GH term. We use the following conven- 
tions: the unit vector n A normal to the brane is taken to point from the region y < 
into the region y > in GN coordinates. The same on both sides of the brane. In terms 
of this vector the induced metric is given by the tangent to the brane components of the 
projection operator gAB = 9AB — n^riB- The extrinsic curvature (the second fundamental 



form of the surface) is defined as Kab = —g^dsDcnD and its trace is K = g Kab- 
Its components in the GN frame are = —^d y g^ u \ y=cons t. For our purposes, we will 
have to evaluate it at y = 0+ and y = 0— , since we put our brane at y = and allow for 
discontinuities of the y derivatives of the metric components. 
The equations of motion for the metric on the brane become 



F(K^ - 9lxv K) + g^F'\ - M 2 FG iiV + M 2 (D^D U F - 9lu/ OF)+ 
X 



+ Xd^d„$ - -g^g a(3 d a ^d^ - g^V = 0. 



(2.10) 



Here □ = g^ u D^D u and is a covariant derivative built upon the induced metric. From 
now on, we shall be using square brackets to denote the discontinuity across the brane 
of the quantity inside, i.e. [W] = W{y = 0+) — W(y = 0—), for any quantity W. Of 
course, if one assumes Zi symmetry across the brane, then W(y = 0+) = —W(y = 0—) 
and [W] = 2W(y = 0+). 

Similarly, the field $ obeys the following equation on the brane 



X' + M 3 F (1) A' 



XW M 2 F^R 

+ XD<D + -—g^d^dv* + 



= 0. 



(2.H) 



3. Static scalar $ = $(y) 



This is the simplest possibility, since to implement the modification proposed in [||, we need 
eventually a non-trivial y-dependent function F. With <I> a function of y only, the system 
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simplifies considerably. Namely, ^-dependent quantities on the brane become constants 
and all derivatives of $ along the brane are zero. In the bulk only y-derivatives survive 
and one is led to the following set of equations. 
In the bulk 



-M 3 FG„ V + M 3 f^F' - \g^gP' - g^P") - y^$' 2 - g^V = 0, 



-M 3 FG yy -^gP' + ^' 2 -V = 0, 



AN' A' _ 
~AN ~ 1 ~ °' 

Where g = g a ^g' a p and □ = d y + \gd y . Contracting the Einstein equations fl3.1|) and 
with g AB one obtains 

3M 3 FR - „ / ~, ~„\ 3X 



2M 3 (gP' + 2F") - ^U' 2 - 5V = 0. 



2 

Under our assumption that P' ^ we can solve (|3.3|) and (|3.5| ) for i? and g, provided 

X ^ — M _ ; . 
^3 F 

This leads to time-independent i? and g. In what follows we assume that condition 
is valid. 

On the brane the field equations are 

-M 3 g^[P'} - M 3 P [K^ - g^K] - M 2 FG^ - g^V = 0, 



[X 1 ] + M 3 PV [K] + M2F(1)R _ V W = o. 



3.1 
3.2 
3.3 

3.4 



3.5 



3.6 



3.7 
3.8 



All tilded quantities are evaluated on the brane and consequently are constants. Contract- 
ing equation (|3.7| ) with g^ v one obtains 

-4M 3 [F'] + 3M 3 P [K] + M 2 FR - W = 0. (3.9) 

This, together with (|3.8| ) form a system of two linear algebraic equations for [K] and R, 
which has a solution provided 3FF^ — 2F^F ^ 0. In this case both [K] and R are 
constants. If, on the other hand, 3FF^ — 2F^F = 0, then for consistency we have to 
require 3F (VW - [X']) = 4FW (v + M 3 [P'}) and [K] and R would not, a priori, have to 



be time-independent. However, given that in our notation K(yo) = —\g^ v <?L, = — \g\ y=yo , 
which is time-independent, we conclude that R is also necessarily time- independent, a 
consequence of our assumption that <I> = &(y). 

It is rather straightforward to take into account the remaining equations. One is led 
to the following conclusions: 
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• R, g, R and [K] do not depend on time 

• A = a(y)a(t) and N = v(y)n(t), i.e. the functions A and N are factorized. 

• There are two possibilities for these functions: either N = b(t)A, or A = and N 
undetermined. 

The second possibility gives Minkowski brane, which is not particularly interesting for our 
purposes. The solution in this case can be found analytically, reduced to quadratures. 
The first possibility, to which we concentrate next, looks more interesting for cosmology, 
because it leads to a time-dependent metric on the brane. 

4. Construction of the solution for static scalar and N = b(t)A 
4.1 Bulk part 

In this case n(t) = b(t)a(t) and a(y) = v(y). The non-zero components of the Einstein 
tensor become 



v" 



G tt = (-b 2 A 2 )3P, G ab = (A 2 5 ab ) (P + 2Q) , G yy = 3 [P + Q - 2— 

1 v' 2 v" n 'b u' 2 v" 

where P = — -^-r H =- H , Q = — ^— H =- H . 

A 2 b 2 v 2 v AAb 3 v 2 v 

Let us, now, take a closer look at the {tt} component of equation (|3.1| ). All terms in this 
equation, with the exception of the first one, are of the form d x (a function of y). This 
is a direct consequence of the factorization property of gtt and the t— independence of <!>. 
Only the first term may a priori depend on time, due to the presence of the —l/(A 2 b 2 ) 
term in P. Thus, the {tt} component of equation Q3.1| ) for non-vanishing d takes the form 
a(t)b(t) = (function of y), and leads to b = l/(Ha) where if is a constant. The Einstein 
tensor then becomes 

G>j, = ^(-H 2 + v' 2 + vv")~g^, G yy = ^(~H 2 + v' 2 ). (4.1) 

The Einstein equations read 

A(-F 2 + v' 2 - vv") = ( X& 2 + M 3 F" - M 3 -F') , (4.2) 

M 3 F V v ) 

: (-H 2 + u' 2 ) = -J-^r (--X& 2 + V + 4M 3 -F') . (4.3) 



6 { _ R 2 ^ = _J_^ + V + 4M 3 ^. 

u 2y M 3 F\ 2 v 

Notice that a(t) remains arbitrary and will not be fixed by the equations of motion. How- 
ever, any form of a(t) leads to (anti) de Sitter brane metric, because it can be absorbed 
into the definition of time by a/{HA)dt = dr, in terms of this new time r one obtains 
a(r) = aoe Hr . So, the bulk metric reads 




~9AB= \ e 2HT / (4.4) 
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where, without loss of generality, ao was chosen equal to 1. 

Given the functions F(y) and v(y) one can use fl4.2|) to obtain X& 2 as a function of y 



X& 2 = - M 3 ( F" - —F' ) (4.5) 



v 

where h = —H 2 + v' 2 — vv" . Finally, plug this into equation ( [4.3| ) and express the potential 
V 

2v 2 2 V v ) v ' 



One can check by direct substitution that equation (2.4) is satisfied. 



4.2 Life on the brane 

The geometry and dynamics on the brane, using the time coordinate r, are given by 

d/iv = \ n oHtz I ' ^ = 12i^ 2 , = -3H 2 g^ u , 
V U e o ab J 

[K^\ = -[v']g^ [K] = -A[ V '\. 

As usual, [z/] = v' + — v'_ is the discontinuity of the function v'{y) across the brane. Z2 
symmetry would force them to satisfy v' + = —v'_. Also we choose i^(0) = 1. 

As discussed above, we have to distinguish two cases, (a) If 3FF^> - 2F^F / wc 
may unambiguously express and H through F, F and V as follows 

F (V« - [X']\ - 2F« (v + M 3 [F'} 



H = ±* 



2M 3 (3FFW - 2FWF) 

(4.7) 



3F (VW - [X'U - 4F(!) (V + M 3 [F'} 
\ 6M 2 (3FFW - 2F( 1 )F) 



(b) If, on the other hand, 3FF^ - 2F^F = 0, we have to require 3F [V^ - [X'\ 
4F« (v + M 3 [F']Y In this case, one of the quantities \u'\ and H remains undetermined. 



However, they are related by 



As a special case, it is easy to check that taking F = F = 1, X = X = V = 0, V = a and 
i/_ = — v'_ = eH one reproduces the known relation ( |1.3| ) between brane-tension and H of 
the DGP setup. 
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5. Application to the modified DGP model and discussion 



Let us recapitulate the steps one has to take to construct a modified model of the type we 
are proposing here, based on <3? which depends only on y. For our purposes, it is convenient 
to start with a given F, since the modified DGP model |l| gives some idea about the 
desirable form of F(y). Namely, F = 1 — m5{y) such that F(y = 0) = const < and 
F(y 3> a) = 1, where a is the width of localization of F near the brane. The metric has 
the general form (|4,4| ). Given, in addition, the function v{y), one uses equations ( |4.5| ) and 
(|4.6| ) to obtain X and V as functions of y. Finally, on the brane, one has to satisfy ( |3.7j ) 
and (|3.§|). It remains to determine the $ dependence of the bulk quantities. For that, one 
needs to specify (or know independently) either <5(y) or the dependence on $ of one of F, 
X or V. 

Consider the special case of Minkowski bulk. This fixes v(y) = 1 + viy, with vi = eH 
and we assume Z 2 symmetry. Take F = 1 — ^^^rr^rj and $(?/) = y- One is led to the 
following model 

~ m 
F = 1 



a cosh (<I> 2 /a 2 ) ' 



mM 3 ( 9/ , , 9, 9nn 9 sinh (2$ 2 /a 2 ) 

X = T — — — 2$ 2 3 -cosh 2$ 2 /a )) +« 



V 



a 5 cosh 3 ($ 2 /a 2 ) y v v ' // 1 + y 

mM 3 / n , , / ^9, 9\\ 9 1 + 8^1$ , / «9 / 9N 

/ ^T%^ /o „„„U /0/T>^ 1 L /O^Tv^ /„.^ 



(5.1) 



— ; 5 „ — 5- ( * (3 - cosh (2$ 2 /a 2 )) + a 2 — ^ i— sinh (2$ 2 /a 

a 5 cosh 3 ($ 2 /a 2 ) V 2(l + i/i$) V 

F = 1, V = a. 



It is easy to check that bulk equations are satisfied. Equation ( |3.8| ) on the brane is satisfied 
trivially, while equation (|3.7|) reduces to 



-QM 4 [l-T.) e H + 3M 2 H 2 - a = (5.2) 



which determines H as 
H 



which coincides with one given in (1.5) up to the rescaling of a by tt and thus explicitly 
realizes the trick outlined in |l]]. 

As another special model, let us take X = 0, still with Minkowski bulk. Then expres- 



sion (4J3) becomes an equation for the function F with the general solution 

F{y) = C l {l + u l yf + C 2 . (5.4) 

Here we see that F is not localized near the brane and becomes important in the whole 
range of y up to the Rindler horizon y = —\jv\. However, the potentially dangerous 
domain of y with F < can be made arbitrarily narrow, with \F(y = 0)| also infinitesimal. 
Indeed, adjust F(y = 0) = C 2 + C\ = —a < to be extremely small. Then, the width of 
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the negative domain becomes — a/(C2^i). In order to have realistic H we have to adjust 
(«M 3 ) ~ (lOOMeV) 3 |§. 

Let us briefly summarize the results of the present analysis. A covariant implemen- 
tation of the modified DGP model M was constructed, by means of a bulk scalar field, 
with appropriate coupling to gravity. Flat metric in the bulk is a possible solution with 
non-constant scalar field. The parameters of the theory can be adjusted in such a way that 
one has not growing metric in the bulk (the usual normal-branch behavior) with positive 
Hubble constant on the brane, leading to accelerating expansion on the brane, even with- 
out any extra matter and even with negative brane tension. We explicitly demonstrated 
two models satisfying our requirements (a)-(d) of the Introduction, (i) In the first one the 
(5-function-like profile of the function F is localized near the brane and in the limit a — > 
the width of localization goes to zero. Thus, in the a — » limit the bulk is modified in 
a narrow domain near the brane only. However, the presence of the kinetic term for the 
scalar field in the bulk is necessary, (ii) In the second model, on the other hand, the scalar 
field does not have a kinetic term, but the profile of F does not go rapidly to unity away 
from the brane. However, one may adjust integration constants and make the y-domain 
with negative F infinitesimal, if necessary to be consistent with observations. 

A detailed stability analysis of the hereby proposed class of models is an open question. 
Also, time-dependent bulk scalar field solutions, which were not studied here, may give rise 
to interesting phenomena on the brane. Finally, asymmetric brane- world setups in the spirit 
of the recent paper [p7g| ] may provide further possibilities for model building. 
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